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Context. The detection of the diffuse supernova neutrino background may be imminent, but theoretical predictions are affected by 
substantial uncertainties. 

Aims. We calculate the signal and its uncertainty with the present configuration of Super-Kamiokande and consider the possibility of 
lowering the threshold by means of gadolinium loading. 

Methods. We model neutrino emission following the analysis of SN1987A by Pagliaroli and collaborators 2009 and use the number 
of expected events in the neutrino detector as a free parameter of the fit. The best-fit value of this parameter and its error are evaluated 
by means of standard maximum likelihood procedures, taking into account properly the correlations. 

Results. The uncertainties in the astrophysics of the emission dominates the total uncertainty in the expected signal rate, which 
conservatively ranges from 0.3 to 0.9 events per year and from 1.1 to 2.9 with gadolinium. 



1. Introduction 

Massive stars end their life exploding as core collapse super- 
novae (SN) and leaving compact remnants, such as neutron stars 
or black holes. Most of the binding energy is released by means 
of neutrino emission, though the details of the emission are not 
yet completely understood. The SN1987A event provided evi- 
dence that the neutrino emission exists and is detectable (Hirata 
et al. [19871 Bionta et al. [19871 Alekseev et al. [19881 . The neu- 
trino signal olfers us the most promising opportunity to date to 
probe this unique astrophysical system. 

Neutrino emission from core-collapse supernovae that have 
exploded in the Universe cumulate in the diffuse supernova 
neutrino background (DSNB). The prospects of detection seem 
promising; for reviews and references to the original papers, we 
refer Ando & Sato (12004) and Beacom (2010l. To predict the 
DSNB signal, two different quantities are needed: the explosion 
rate of core collapse SNe as a function of the redshift and the 
average neutrino emission of the individual supernovae. Each of 
these quantities is known up to an error which implies uncer- 
tainty in the prediction. Until both uncertainties are quantified, it 
is impossible to evaluate the reliability of the expectations. 

The rate of core-collapse supernovae can be identified with 
the rate of massive-star formation as a function of redshift. This 
quantity, in turn, can be obtained from knowledge of the cos- 
mic history of star formation and the initial distribution of mass. 
Horiuchi et al. ( I2009I I used a comprehensive compilation of data 
to evaluate the uncertainty in the rate. 

We attempt is to complement the study of the uncertainties 
by evaluating the impact of the neutrino emission for an individ- 
ual supernova. Following Fukugita & Kawasaki (I2003I I. we base 
our inferences on SN1987A observations using for the analy- 
sis the neutrino emission model of Pagliaroli et al. (2009). We 
evaluate the expected signal and its uncertainty for the Super- 
Kamiokande detector. We consider its present configuration, 
22.5 kton of fiducial volume with a threshold of E,/, = 19.3 MeV, 



Malek et al. (120031 ) and consider the possibility of lowering the 
threshold to 1 1.3 MeV by employing gadolinium loading as ad- 
vocated by Beacom & Vagins (i2004i) . 

2. The emission from an individual supernova 

We determine the spectrum from SN1987A data, which we 
then use for the analysis of DSNB. This is done by parameteriz- 
ing the time and energy distribution, then by fitting the data from 
Kamiokande-11, 1MB, and Baksan and finally by integrating over 
the time distribution. The main motivation behind this procedure 
is that we still do not have a definitive theory of supernova ex- 
plosion; thus, despite their paucity, the data from SN1987A play 
a very important role in guiding our understanding of supernova 
neutrino manifestations. 



The model for supernova emission Motivated by the prospect 
of exploring, through the neutrino signal, the physics and astro- 
physics of the gravitational collapse, we proposed in Pagliaroli 
et al. ( 120091 ) a parameterization of Vg emission, based on our 
present understanding of emission processes improving the 
model by Loredo & Lamb (i2001 ). Our model has two emission 
phases and the Vg flux is 



<l),,(f, £') = (5,(f, £') + [!- jk(t)] ^c(t - T„, E'), 



(1) 



where t is the emission time and E' is the emitted neutrino en- 
ergy. The first term is the flux generated during the phase of 
accretion and above the shock by the interactions between the 
positrons and the target neutrons. It describes volume emission. 
It uses three free parameters: the initial accreting mass M^, the 
time scale r,,, and the initial temperature of the positrons Tg. The 
second term <^c is the flux coming from the thermal emission of 
the new-born proto-neutron star (i.e., the cooling phase). It de- 
scribes radiation from a surface, whose area is proportional to 
the radius of the neutrino sphere Re, with a time scale t^. and the 
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Fig. 1. The SN1987A best- fit spectrum (continuous line) and the 
two contributions from accretion phase (dotted Hne) and cooling 
phase (dashed line). In the inset, we plot the percentage differ- 
ence with the Fermi-Dirac approximation of the spectrum, given 
in the text. 




Fig. 2. The cosmological factor n of Eq. Q for various fits of the 
core-collapse rate. The continuous line is for the fiducial rate, 
the dotted line is for the upper one, and the dashed line is for the 
lower rate. 



initial temperature of the emitted antineutrinos Tc. Finally, the 
function jkit) links smoothly the two emission phases, delaying 
the cooling emission by Tg. Analytic expressions of these three 
functions, Oa, O^, and jk{t) are given in Eqs. 10, 13 and 18 of 
Pagliaroli et al. (120091) . respectively. We include neutrino oscil- 
lations with normal hierarchy as discussed in Sec. C of Pagliaroli 
et al. (|2009] l. 

SN1987A data analysis In Pagharoh et al. ( I2009I I, we tested 
the parameterized model using the small set of events collected 
in 1987 by Kamiokande-11, 1MB, and Baksan, leaving aside 
an interpretation of the events recorded by LSD (Aglietta et 
al. I1987i When we compared the fit based only on the cool- 
ing phase, the one adopted in usual SN1987A data analyses, 
with our fit, we found that the two-phase emission model is 50 
times more probable (i.e., we got a 2.5cr indication in its fa- 
vor). Moreover, the best-fit values of the astrophysical param- 
eters, namely Rc = 16 km, M„ = 0.22Mo, = 4.6 MeV, 
Ta - 2.4 MeV, Tc - 4.7 s, and Tg - 0.55 s, agree well with 
the general expectations. For instance, the duration of the accre- 
tion phase is shorter than one second, the radius of the neutrino 
sphere is similar in size to the neutron star, and the total radiated 
energy 2.2 x 10^^ erg is similar to the binding energy. Finally, the 
luminosity curve and the mean energy as functions of time both 
resemble the results of numerical simulations. The errors in the 
parameters and their correlations are large, as expected from the 
limited number of events (see again Pagliaroli et al. 2009). 

The previous analyses of SN1987A data of Lunardini (2006i 
and Yuksel and Beacom (i2007i) designed to predict DSNB fo- 
cused only on the energy spectrum, motivated by the opinion 
that the time distribution of the events is irrelevant. However, we 
show that our detailed theoretical description of neutrino emis- 
sion leads to a peculiar integrated spectrum that retains an im- 
print of the two different emission phases. 

The emission spectrum By substituting in the Ve flux of Eq.® 
the best-fit values reported in the previous paragraph and inte- 
grating both in a 30 s time window and on the emitting surface, 
we obtain the reference spectrum dN/dE' for a single supernova 
event from SN1987A shown in Fig. [1] We compare this spec- 
trum with a Fermi-Dirac distribution, i.e., 4r-r. — tpttt — To ob- 
tain the same total integral and the same first two momenta be- 
tween our spectrum and the Fermi-Dirac distribution, we need 



k - 2.97 X 10 , a temperature T — 3.76MeV, and a pinching 
factor J] = 0.679. However, the percentage difference between 
the Fermi-Dirac spectrum and our spectrum is not negligible as 
shown in the inset of Fig. [T] The agreement worsens using the 
parametrization of Keil et al. (2003). The emitted electron an- 
tineutrinos of our best-fit spectrum are of relatively low energy; 
half of them are emitted at energies below 11.2 MeV. In particu- 
lar, when we calculate the number of events A^^v expected from a 
galactic supernova assuming an energy threshold of E,„i„ > 6.5 
MeV,A^,, oc JdE'cr{E)<:>y^ as a function of the upper extreme of 
integration E,„ax, we find that there are the 25%, 50%, and 75% 
of the events below E,„ax-14; 18, and 24 MeV, respectively^ 

3. Expectations for the diffuse neutrino flux 

We estimate the diffuse neutrino flux accumulated by all the past 
supernovae exploded in the Universe, assuming that the antineu- 
trino flux discussed in the previous Section represents the typical 
emission of a core-collapse event. 

In a Friedmann-Robertson- Walker flat Universe, the ex- 
pected DSNB flux is 

dE Ho Jo ^ ^JQ.J^ + Q,„(l + zy ^E' 

The last term of the integrand, dN/dE', is the spectrum of a 
single SN emission discussed above and calculated for the red- 
shifted energy E' - {I + z)E. The Hubble constant is Ho = 
71 km s^'Mpc"', c is the light speed, and the values of the cos- 
mological constants are fl,„ - 0.27, D.^ - 0.73 as measured 
from the WMAP experiment in Jarosik et al. ( 1201 Oi l. These as- 
sumptions do not introduce significant errors into the predic- 
tions. 

The key term we have to discuss is RccsNiz), namely the 
rate of core-collapse SN for comoving volume as a function of 
the redshift z. This function can be obtained from the number 
of stars formed in a comoving volume, i.e., from the star for- 
mation rate as a function of the redshift, and from the fraction 
of stars in the appropriate mass range for core collapse SNe. 

' Of course, cr denotes the cross-section for tiie IBD process v^p — » 
e*n. We use Eq. 25 of Strumia & Vissani ( I2003t and tiirougiiout tiie 
paper, we adopt tiie approximation Ey^ = Eg* + 1.3 MeV, wliicli is ad- 
equate for our needs. Unless specified otherwise, we always refer to 
energy. 
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Fig. 3. Expected DSNB flux in a logarithmic scale for the three 
different fits of core-collapse rate. Also is shown the upper limit 
of this quantity given by Super-Kamiokande in 2003 in the en- 
ergy range between 19.3 and 27.3 MeV. The meaning of contin- 
uous, dotted and dashed lines as in Fig. (|2]l. 



Following Horiuchi et al. ( 120091 1. we adopt the initial mass func- 
tion of Salpeter (1955') and three different analytic fits for the 
rate of core-collapse to take into account the astrophysical un- 
certainty. In Fig.©, we show the quantity 



n(z) = 



VnA+0„,(l+Zp 



(3) 



which we call the "cosmological factor", namely the first in- 
tegrand term in Eq. This plot summarizes all the relevant 
assumptions about the cosmology and the distribution of the 
sources. Horiuchi et al. (i2009) demonstrated that the use of ei- 
tiier tiie Kroupa (l200TT l or Baldry & Glazebrook (12003^ initial 
mass functions, or the uncertainties in the lowest mass that forms 
a supernova, introduces only a negligible error in the estimates. 

We now use Eq. (|2]i to calculate the flux of the diffuse SN 
neutrino background d<p/dE as a function of the neutrino energy 
at the Earth. We plot it in Fig. (O along with the upper bound 
in the energy region E > E,!, = 19.3 MeV of 1.2 Vp cm"^ s"' 
obtained by the Super-Kamiokande collaboration (Malek l2003l l. 
The total flux of DSNB with the SN1987A best-fit model is = 
27.2 cm"^ s"' for the fiducial rate. Only ~ 1 % of this flux can be 
actually observed with a neutrino energy threshold of 19.3 MeV; 
with the use of gadolinium, the accessible spectrum becomes 
~ 8%, showing that most of the flux falls in the very low energy 
region. 

Event rate in Super-Kamiokande We consider the events ex- 
pected in a detector such as Super-Kamiokande, with a fidu- 
cial mass of Md = 22.5 kton of water and a detection effi- 
ciency e set equal to the 98% above a neutrino energy threshold 
(Eth = 19.3 MeV in the present configuration or £,/, = 11.3 MeV 
by loading the detector with gadolium). However, we note that a 
scintillator-based detector with average chemical formula C9H21 
(or C6H3(CH3)3) and with 16 (or 25) kton of mass has the same 
number of target protons A^^ = 1 .5 x lO"'-', thus the same number 
of Vep — > e^n interactions. The event rate is calculated easily as 



,.^Np f dE(T(E)^6(E). 



(4) 



For the best-fit model adopted, we expect to see in Super- 
Kamiokande the number of events per year reported in TabledTJ 
for two different energy thresholds and three different fits of the 
rate. We expect 0.39 - 0.65 events per year when we consider 
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Fig. 4. Ax^ curves as a function of event rate Nev in Super- 
Kamiokande, as obtained from SN1987A data analysis and for 
two different energy thresholds. The set of thin lines correspond 
to the present threshold of 19.3 MeV whereas the set with thick 
lines correspond to £,/, = 11.3 MeV. The meaning of the indi- 
vidual lines is as in the previous figures. 



a threshold of 19.3 MeV; this increases to 1.35 - 2.35 events 
per year lowering the energy threshold to 11.3 MeV. Our results 
agree well with those of Horiuchi et al. ( I2009I I. 



4. Uncertainty due to the emission model 

In Sect. [3j we discussed the expectations for the DSNB flux 
based on the SN1987A best-fit model and we considered only 
the uncertainty related to the core-collapse rate. However, the 
SN1987A data set is small and the procedure of using only the 
best-fit model is doubtful, until one estimates the associated the- 
oretical error. The main aim of this work is, indeed, to quantify 
this uncertainty, which is an entirely new result. 

We proceed as follows. The likelihood in the SN1987A data 
analysis used by Pagliaroli et al. (2009') is a function of six un- 
known astrophysical parameters. The same is true for the DSNB 
rate A^^.,,. For a given energy threshold, this can be written ex- 
plicitly as Nev = Ma fija, Tg) + gijc, Tc, T^), whcrc the two 
contributions correspond to the two terms of Eq. ([1] and / and 
g are known functions. The idea is simply to rewrite the likeli- 
hood, substituting the neutrino-sphere radius Rc with the quan- 
tity in which we are interested, A^^,,,. We can then apply a standard 
marginalization procedure: for each fixed value of N^v and vary- 
ing the other five parameters, we find the maximum value of the 
likelihood L. Thus, we obtain Ax^ = 2(ln L,„a.x - In L) a function 
of only. 

The results of the numerical evaluations of this Ax^ are 
shown in Fig.©. The thin lines represent the neutrino energy 
threshold used in Malek et al. (120031 1. £,/, - 19.3MeV, whereas 



Table 1. Rate of DSNB events expected in Super-Kamiokande 
detector [yr"']. Together with the best-fit values, we show the 
Icr and 2cr statistical errors obtained from marginalization pro- 
cedure of Sect.m 



CCSN Rate 


£•„, = 19.3MeV 


= 11.3MeV 


Upper 
Fiducial 
Lower 


o-5i!!;:l^(i^)!!!1^(2^) 


2.35!«;^»(1<T) *J,jt(2a) 
1.82!«;^«(l(r) !«;«](2tr) 
1.35!«i^(lfr) tJ^(2tr) 
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Fig. 5. Cumulative distribution of the events rate in Super- 
Kamiokande, with the same notation as Fig. (|2]) for the lines. 
We give the 2cr error bars on that prediction when - 11.3 
MeV and = 19.3 MeV. 



the thick ones correspond to the threshold of Eth - 11.3MeV 
that can be achieved with gadolinium. Three hnes of each type 
are given, for the three different fits of the cosmic rate for core- 
collapse events. The Icr and 2cr errors found from the Ax^ curves 
are reported in Table ([1). 

We note that this procedure properly takes into account the 
correlations among the parameters, that are dictated by the like- 
lihood function, without the need to explicitly calculate the co- 
variance matrix. 



5. Summary and discussion 

We have evaluated the signal expected for diffuse supernova neu- 
trino background in Super-Kamiokande and its associated uncer- 
tainties, using SN1987A data to constrain the model of neutrino 
emission. The results are summarized in Fig.(|5), where we show 
the cumulative rate of DSNB events in Super-Kamiokande as a 
function of the energy threshold for the three different fits of the 
core-collapse rate. 

Following Horiuchi et al. ( 120091 ), we used the three different 
fits as a "generous" assessment of the cosmological uncertainty, 
although we are unable to attach to this range a precise statisti- 
cal meaning. The result is given in Table ([TJ: the cosmological 
uncertainty is ~ 25% for the higher threshold and ~ 27% for the 
lower threshold, in reasonable agreement with Horiuchi, despite 
the different emission model used. 

On the other hand, selecting the fiducial cosmological model 
for the rate, we evaluated the uncertainty on the DSNB due to the 
emission model. The Icr (resp. 2cr) percentage error is ~ 33% 
(resp. 65%) for the higher threshold and becomes ~ 23% (resp. 
46%) for the lower threshold. The prediction is more precise in 
the second case since the SN1987A data provide a tighter con- 
straint of the low energy region of the spectrum. On the basis of 
the results of Pagliaroli et al. (I2009I I. we conservatively estimate 
a 10% systematic uncertainty in our spectrum, which is intrinsic 
to the parameterization used in the fit. Our 2cr ranges are curi- 
ously similar to those obtained by Ando & Sato (2004), who use 
three different emission models obtained from numerical simu- 
lations. 

It is difficult to combine the two errors in a robust way, since 
only the one related to the emission model has a precise statisti- 
cal meaning. We thus construct a global range for the expected 
rate of DSNB events by the following conservative procedure. 
The upper value of the range is obtained by summing the best- 
fit value for the upper core collapse SN rate and its upward 1 cr 



(resp., 2cr) statistical error; similarly, for the lower value of the 
range. In particular for £,/, - 19.3 MeV, the expected event rate 
in Super-Kamiokande ranges between 0.27 (resp., 0.19) events 
per year and 0.88 (resp., 1.16) events per year, giving a global 
uncertainty of ~ 53% (resp., 95%). For the lower energy thresh- 
old Eth - 11-3 MeV, the range becomes 1.06 - 2.94 (resp., 
0.83 - 3.59) events per year with a percentage total error of 47% 
(resp., 76%). As quantified by the percentage errors, most of the 
global uncertainty is due to the emission model uncertainty and 
to the low number of SN1987A data. 



Appendix A: An alternative expression for the rate 

We introduce an alternative expression for A^^,, that emphasizes 
the role of the supernova flux at the source providing additional 
insigh into the results. From the emission spectrum dN/dE', we 
define an effective flux 



dE And^T dE' 



with ['-4^-^'^^^^ 



(A.l) 



where we have introduced the typical value for the cosmic den- 
sity of «, - 2x 10"^ SN/(Mpc-' yr) and the observational time 
T . By using as an integration variable E' - E{\ + z) (i.e., the 
neutrino energy at the emission), we rewrite the signal rate in 
Eq. dUi mimicking closely the signal from a galactic supernova 



(A.2) 



Ne^^Np dE'cr{E')—(E')e,{E'). 
Je,„ dE' 

The distribution of cosmic supemovae is contained in the func- 
tion 

r^'^dz n(z)o-{E'Ki+z)) 

e.(£ )= I — — e(£7(i+z)). (A.3) 

Jo i+z n. o-{E') 

This function leads to a severe cut in the rate of events given 
by Eq. (IA.2b . Above the energy threshold, this function can be 
approximated linearly as 



e*{E') cc E' - Ell,, 



(A.4) 



which corresponds to DSNB selecting the highest energy tail of 
the spectrum. This is important in connection with SN1987A. 
We have only limited information on the highest energy tail of 
the spectrum, mostly thanks to 1MB, as noted by Fukugita et al. 
2003 We can constrain far more reliably the spectrum at lower 
energies, thanks to data from Kamiokande-11 and Baksan. Thus, 
we expect, by lowering the threshold, not only the number of 
events due to DSNB to increase, but also the accuracy of the 
predictions based in SN1987A to improve, in agreement with 
Fig. |5] We also note that the numerical value of d in Eq. (lA.lb 
agrees with the DSNB signal being modest. 
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